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Given O/n in Aj the routine places (l/2} sin in A 
when a standard entry is usedo 

(1) If A contains l/2 + ©A* the qmuitity 1/2 cos d 
will be obtained. 

(2) If d/jc overflows in A from left shifts or addition 
the result is unaffected because sin (O 4- 2n it) = sin 0« 
(l) If memory location q contains 9/^, (l/Z) cos 9 may 
be obtained vith the entry 

LJ q 
50_£ 



P 

p + 1 

(2) To find sin 3*32 radians when memory location q 
contains O.3B2 and memory location r contains lO/^^ we 
proceed as fo3JLow8, ignoring the overflow in A» 
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(1) Given Q in radians we evaluate sin (0/2) by a 
Tchebyscheff polynomial expansion to the Hfeylor series 
expGmsion of the sine fimction^ using Libraiy routine 
KA-1. 



(2) If 9 is scaled so that -1 < 2""^9 < 1^ we form 
2^ (2-^0/jt = 0/rt (mod 2.)) = 0»/trtc 

(3) Using the identity sin a = sin (ic - a) we obtain 
9'* with °l/2 < 9"^ < 1/2 such that sin*«" = sin 9' 
= sin 9» 

ERROR ANALYSIS (l) In word k we develop -k (q/ti)^ , Consider 

|9/j(| < 1/2. This quantity, possesses up to 58 significant binary digits. We 

form (o/jt) < 1/^ and shift left twice^ introducing thereby at most an error 
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of 2 _, since q contains a significant binary digit, 

(t) We next develop the polynomial: 

1/2 sin 9 = Q/yt [a (29/ir)-^ +^a^^(29/7i)^^ + 

aQ(20/:t)^ +r.a„(29/jt)^ + a^(29/jt)^ + 



a^(29/jtf + a^] 
This is accomplished by the recursion relation: 



2 (1=0,1,2,... 6) 

At each one of these steps we introduce a round°off 
error due to the multiplication. Letting (29/it) = x^ we obtain: 

-±72- sm 9 = — La^ ^x + a, ^x + a^ + a^ ^ a^x + a^x + a, 
jtl3 119'(5 31 

+ 2-^°(x5 + x\+ x5 + x2 + X + 1) eA + 2-^°J 

Since x cannot exceed 1 and |9/it| < l/S we have as error 

-,1^^ -IfO -=.38 

3.2 +2 • The maximum total generated error is therefore 2 ^^'^ 

(3) There will also be developed a propagated error due 
P 

to the inaccuracy of the quantity 4(9/it) . We will admit only the error result- 

P 
ing from the term a-,(29/jt) in equation (l) due to the relative small contribu-= 

tions propagated by the other terms o 

_ 2 - 



Thias: 



a^(20Af = M^ ^1-^ + 2"^^] = .4(29// + (.1^)2"^. 



The propagated error is therefore s (.4)2" ''^ 

(4) The total loaximum error is therefore s 
2^58 ^ (.i^)2"59 = 2 + .4) 2-59 ^ (2,i,)2-'59 

(5) A sample of calculations for the sine of a nxmiber 

-12 
of arguments was found to be in error at most 5 x 10 

which is less than the maximum error indicated above. 

(6) Finally^ we consider the case [Q/^tl = I/2. The 



error was found experimentally to be less than 2 
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Prepare round off 
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1/4 sin 9 
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1/2 sin 9 
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-a^^a = 01797591 X io"5 
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a^2 = *80219199 X 10"^ 
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_a„/2 = .2540876^^99 X 10"^ 
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at/2 = .398ii-63l500ii- x lO"-*- 
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a^/2 = .7853981633^7 
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